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ABSTRACT 



We present multi-configuration Breit-Pauli AUTOSTRUCTURE calculations of distorted-wave photoionization (PI) cross sections, 
and total and partial final-state resolved radiative recombination (RR) and dielectronic recombination (DR) rate coefficients for the 
first six ions of the trans-iron element Se. These calculations were motivated by the recent detection of Se emission lines in a large 
number of planetary nebulae. Se is a potentially useful tracer of neutron-capture nucleosynthesis, but accurate determinations of its 
abundance in photoionized nebulae have been hindered by the lack of atomic data governing its ionization balance. Our calculations 
were carried out in intermediate coupling with semi-relativistic radial wavefunctions. PI and recombination data were determined 
for levels within the ground configuration of each ion, and experimental PI cross-section measurements were used to benchmark our 
results. For DR, we allowed An = core excitations, which are important at photoionized plasma temperatures. We find that DR 
is the dominant recombination process for each of these Se ions at temperatures representative of photoionized nebulae (~10'' K). 
In order to estimate the uncertainties of these data, we compared results from three different configuration-interaction expansions 
for each ion, and also tested the sensitivity of the results to the radial scaling factors in the structure calculations. We find that the 
internal uncertainties are typically 30-50% for the direct PI cross sections and ~10% for the computed RR rate coefficients, while 
those for low-temperature DR can be considerably larger (from 15-30% up to two orders of magnitude) due to the unknown energies 
of near-threshold autoionization resonances. These data are available at the CDS, and fitting coefficients to the total RR and DR 
rate coefficients are presented. The results are suitable for incorporation into photoionization codes used to numerically simulate 
astrophysical nebulae, and will enable robust determinations of nebular Se abundances. 

Key words, atomic data - atomic processes - HII regions - nucleosynthesis, abundances - planetary nebulae: general - stars: evolution 



1. Introduction 

The abundances of neutron(n)-capture elements (atomic number 
Z > 30) are sensitive probes of the nucleosynthetic histories of 
stellar populations and the chemical evolution of the Univers e 
(IWallerstein et al.l 'l997: Busso et al. 199^ ISneden et al.ll2008l) . 
Our knowledge of the production of these elements is largely 
based on the interp retation of abundances derived from stellar 
spectroscopy (e.g.. 'Smith & Lambert" 1990': 'Busso et al.| [T999 
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[Xravaglio et al. 2004; Snedenetal. 2008; Roederer et al. 
and references therein). However, only a limited number of 
trans-iron elements can be detected in stellar spectra. Moreover, 
some classes of stars or stages of evolution with high mass-loss 
rates are very difficult to study spectroscopically due to photo- 
spheric obscuration. 

The recent detection of emission lines from n-capture 
element s in a la rge nu mber of planetary nebulae (PNe) 
(Sterlin g etal .11200 7; Sharp ee et cd]|2007l: [Sterling & Diners teini 
r2008tlSterUngeta l. 2009), as well as other astr ophysical objects 
inclu ding H II regioi is (e. g., Ascin 1994; Lumsden & Piixleyl 
1996 ; Luhman et al] 119981; Baldwin et al. 2000; Puxley et al. 
2000i; Okumuraet al.' '200 1 ; Blum & M cGregor. i2008t 
Roman-L opes et al. 2009), starburst galax ies (IVanzi et al.l 
2008 *). and active galactic nuclei (Thomp son et al J Il978h . 
has demonstrated that nebular spectroscopy is a potentially 
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powerful tool for investigating the nucleosynthesis and chemical 
evolution of trans-iron species. The abundances of «-capture 
elements are of particular interest in PNe, since these species 
can be produced in PN progenitor stars via slow n-capture 
nucleosynthesis (the ".s-process") during the asymptotic giant 
branch (AGB) stage of evolution. 

Nebular spectroscopy provides access to many elements that 
cannot be detected in cool giant stars, such as the lightest n- 
capture elements (Z = 31-36) and noble gases. Due to the dif- 
ficulty in detectin g these sp ecies in stellar spectra or supernova 
remnants (Wallerstein et al.iri995 ). the origins of light n-capture 
elements are based predominantly on theoretical considerations 
that lack empirical validation. In addition, nebular spectroscopy 
enables investigations of nucleosynthesis in classes or evolution- 
ary stages of stars obscured by optically thick circumstellar en- 
velopes. For example, intermediate-mass AGB stars (4-8 Mq) 
experience heavy mass loss, shielding their photospheres in a co- 
coon of circumstellar material that hinders spectroscopic investi- 
gations in the wavelength regimes of many useful n-capture ele- 
ment transitions. As a result, the contribution of these stars to the 
Galac tic inventory of heav y element nuclei is poorly understood 
(e.g.. iKarakas et alj |2009|) . However, these intermediate-mass 
stars produce PNe in which n-capture element emission lines 
are readily detected (Sharpee et al. 2007; Sterling & Diners teinI 
2008). Likewise, AGB stars that become carbon-ric h as a re- 
sult of convective dredge-up (see iBusso et al ] ll999h are char- 
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acterized by high opacities (lMarigoll2002l) that substantially in- 
crease their mass-loss rates, and hence i-process enrichments 
engendered by the final stages of AGB evolution are not well- 
constrained. PNe are composed of material from the stellar en- 
velope at the end of the AGB phase, and hence are useful probes 
of ^-process enrichments during late AGB evolution. 

However, the accuracy of nebular n-capture element abun- 
dances is hindered by the poorly known atomic data for these 
species, thereby limiting the usefulness of nebular spectroscopy 
for studying the origins of these elements. The reason for this is 
that generally only one or two ions of n-capture elements have 
been detected in individual nebulae, and hence the abundances of 
unobserved ions must be estimated in order to derive elemental 
abundances. This is most robustly achieved by numerically sim- 
ulating the thermal and ionization str ucture of nebulae w ith pho- 
toionization codes such as Cloudy ( Ferland et al.lll998h . but the 
reliability of these models strongly depends on the availability 
of accurate atomic data for processes that control the ionization 
balance. For photoionized nebulae such as PNe, these data in- 
clude photoionization (PI) cross sections and rate coefficients for 
radiative recombination (RR), dielectronic recombination (DR), 
and charge transfer (CT). Unfortunately, such data have not been 
determined for the vast majority of trans-iron element ions. 

To address this need, we have computed multi-configuration 
Breit-Pauli (MCBP) distorted-wave PI cross sections and RR 
and DR rate coefficients for the first six S e ions, using 
the atomic structure code AUTOSTRUCTURE ("Badneir'igS^ 
Ll997). Along with Kr and Xe, Se is one of the most widely 
observed w-capture elements in io nized nebulae jSharpee et al.l 
120071; ISterling & Dinersteinll2008l) . In fact, Se has been identi- 
fied in nearly twice as many PNe as any other trans-iron element, 
and hence is the initial target for our study. In subsequent papers, 
we will present similar data for low-charge Kr and Xe ions, as 
well as CT rate coefficients for several «-capture elements. These 
data will be suitable for incorporation into photoionization codes 
such as Cloudy, which can be used to robustly determine cor- 
rection factors for the abundances of unobserved ions, enabling 
much more accurate abundance determinations of «-capture ele- 
ments in ionized nebulae than previously possible. 

Until our study, the photoionization and recombination prop- 
erties of Se ions had received very little attention, aside from 
a handful of PI s tudies of neutral Se (|Manson et aTl 1 19791: 
iGibson et al. 1986; Ch en & Robicheau x 1994*). which is a trace 
species in ionized nebulae such as PNe. We therefore utilize re- 
cently co nducted experimental PI cross-section measurements of 
Se ions (Esteves et all 120091; ISterling etaO 120111; lEsteves et all 
12010 , 201 1) as a comparison and benchmark to our calculations. 

We estimate the uncertainties in our computed data by utiliz- 
ing three different configuration-interaction (CI) expansions for 
each ion, and test the sensitivity of our results to various internal 
parameters of the AUTOSTRUCTURE code. Upon completion 
of our theoretical study of Se, Kr, and Xe ions we will use Monte 
Carlo simulations with photoionization codes to reveal the ef- 
fects of these atomic data uncertainties on nebular abundance 
determinations, illustrating which systems and atomic processes 
require further theoretical and/or experimental analysis. 

This paper is organized as follows; in Sect. |2] we provide de- 
tails of our calculations, and in Sect. [3] we present our computed 
PI cross sections and RR and DR rate coefficients along with 
their associated uncertainties. Finally, in Sect. |4] we summarize 
our results and provide concluding remarks. 



2. Calculations and methodology 

We have computed the electronic structure and MCBP 
distorted-wave PI cross sections for Se°-Se^^ with the 
AUTOSTRUCTURE code (Badnell 1986, 1997). Higher charge 
states were not considered since they are negligibly populated 
in PNe, who se central stars are not sufficiently hot (ge nerally 
< 2 X 10^ K; lNapiwotzkilll999l;IStanghellini et alj|2002t) to sig- 
nificantly ionize species with ionization thresholds greater than 
100 eV. RR and DR rate coefficients were determined from the 
direct and resonant portions of the PI cross sections, respectively, 
using detailed bal ance. We used th e independent processes ap- 
proximation (Pind zola et al.l [T992h to treat RR and DR sepa- 
rately. 

The theoretical background of such calculations ha ve been 
provided elsewhere (iBadnell et al.l2003l:lBadnelll2006bl) . and we 
refer the reader to those references for a full discussion. In each 
of the following subsections, we discuss details of calculations 
for each process considered. 

2.1. Electronic Structure 

To test the sensitivity of our results to the adopted CI expan- 
sions, we constructed three different CI expansions for each Se 
ion (designated "small," "medium," and "large"). The medium 
configuration sets provide the best compromise between accu- 
racy and computational expense, and are adopted as the basis 
for our PI, RR, and DR results. The CI expansions for each ion 
are listed in Table [10 

In all cases, the electronic structure was computed 
with Thomas-Fermi-Dirac-Amaldi model potentials, adopting 
intermediate-coupling an d /c-averaged relativistic wavefunctions 
(ICowan & GriffinI Il976h . We Schmidt orthogonalized the k- 
averaged orbitals in our calculations, but compared our PI and 
RR results with those obtained without forcing the radial orbitals 
to be orthogonal. The average of the LS term energies were op- 
timized by varying the orbital radial scaling parameters so as 
to best reproduce experimental energy levels and ionization po- 
tentials from NIST (Ralchenko et al. 2008). The radial scaling 
parameters adopted for each CI expansion of each ion are given 
in Table |2] A comparison of selected calculated and experimen- 
tal energies is shown in Table [3] We also compare our Einstein 
A-coefficients with those available in the literature in Table |4] 

The structure computed with the medium configuration sets 
generally reproduce experimental energies and ionization poten- 
tials to a good degree of accuracy, lending credence to our use 
of orthogonalized radial orbitals. The ionization potentials are 
within 2-3% of the NIST values for all ions with the exception 
of Se'^, which is a very challenging system to model but fortu- 
nately is negligibly populated in ionized nebulae such as PNe. 
The energies of most levels tabulated by NIST are reproduced to 
within 5% for Se^-Se^+, and to within 9% for the neutral case. 
A few exceptions occur, particularly for levels within the ground 
configuration, where the discrepancies with experiment are as 
large as 23% for Se+ 4p^ ^Dj,/2. 

The accuracies of the level energies aff'ects those of the 
Einstein A-coefficients A,y we computed. For forbidden transi- 
tions within the ground configurations, the main comparisons are 
the theoretical studies of Biemont & Hansen in the 1980s (see 
Table |4] for references), who utilized the Hartree-Fock with rel- 
ativistic corrections (HER) and Hartee-plus-statistical-exchange 
(HXR) methods. In many cases, our computed agree with 



' Tables 1-4 are available in the online version of this article. 
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those previous studies to within 30-50%, but exceptions occur, 
particularly for the lowest-charge states. However, it is worth 
noting that our CI expansions generally are larger than those 
used by Biemont & Hansen, and hence include configuration 
mixing not incorporated in their calculations. 

2.2. Photoionization 

MCBP distorted-wave PI cross sections were computed in inter- 
mediate coupling for neutral Se and each of the first five ions. 
Cross sections were calculated for each level in the ground con- 
figurations, as the vast majority of ions in photoionized plas- 
mas reside in the ground configuration. The CI expansions listed 
in Table [T] were used in these calculations, with radial scaling 
parameters from Table |2] In the cases of Se+, Se^^, Se-'^, and 
Se^^, experimental PI cross sections measured at the Advanced 
Light Source (ALS) synchrotron radiation faciUty at Lawrence 
Berkeley National Laboratory in California were used to test the 
accuracy of our computations (see Sect. 13.1.2b . 

PI cross sections were computed up to 100 Ryd for 
each ion. The length gauge was used at low energies, and 
AUTOSTRUCTURE automatically switched to velocity gauge 
at higher energies. Because of our use of /f-averaged relativis- 
tic orbitals, it was not possible to use acceleration gauge. The 
electron energy mesh was interpolated onto the photon energy 
mesh using a 4-point Lagrange interpolation. Photoexcitation- 
autoionization resonances were included for each system, over 
a more restricted energy range due to the finer energy mesh 
required for comparison to experimental measurements. In all 
cases, we utilized the radial scaling parameters associated with 
the (A^ + l)-electron ion, as tests showed that using the target 
scaling parameters led to very poor agreement with experimen- 
tal PI cross sections. Scaling parameters for continuum orbitals 
were taken to be equal to those of the highest principal quantum 
number bound orbitals with the same orbital angular momentum 
for s, p, and d orbitals, and equal to unity otherwise. 

2.3. Radiative recombination 

Total and partial final-state resolved radiative recombination rate 
coefficients were computed from the direct PI cross sections us- 
ing detailed balance. The CI expansions of the (A^ + l)-electron 
ions were augmented to include one-electron additions to the 
target's ground and mixing configurations. Since it was not pos- 
sible to use acceleration gauge at high energies, the rate coeffi- 
cients may be inaccurate at the highest temperatures. However, 
we have tested that this does not affect the accuracy at photoion- 
ized plasma temperatures (see Sect. 13.2.1b . The rate coefficients 
were determined over the temperature range (10' - 10^ K, 
where z is the charge. We used a fixed maximum principal quan- 
tum number of 1 000, and a maximum orbital angular momentum 
of 4 for these calculations. For / > 4, the analytic hydrogenic ra- 
dial integrals from the recurrence relations of Burgess. (.1964.) 
were employed. 

2.4. Dieiectronic recombination 

In photoionized nebulae such as PNe, low-temperature DR is the 
dominant form of recombination for many species. We therefore 
computed DR for An = core excitations, neglecting An > 1 
excitations, which are usually negligible at the temperatures of 
interest (^lO"* K). However, for DR of Se^ and Se^^ (forming 
Se" and Se^, respectively), we allowed core excitations into the 



5i and 5p orbitals since the lowest 4s^4p*5i and As^Ap^5p 
(k = 2,3) levels lie below the lowest As^ 4/?* 4-d levels in energy 
for these ions. Due to the intensive computational requirements 
of DR calculations, we ignored the (small) contributions to the 
Se^ and Se^^ DR rate coefficients from excitations into the 5d, 
5/, and 5g orbitals. 

In the case of low-temperature DR, it is critical to have ac- 
curate energies for low-energy autoionizing states, in order to 
determine whether they lie just above the ionization threshold 
and contribute to DR, or just below threshold in which case 
they do not. Unfortunately, the energies of these levels have not 
been spectroscopically determined for Se or other n-capture el- 
ements, thereby limiting the accuracy of our calculated DR rate 
coefficients. To alleviate this source of uncertai nty, we used ex - 
perimental target energies whenever possible ( Badnellll2006al) . 
When experimental energies are incomplete for a given term, we 
assumed the theoretical level splitting. We did not attempt to ad- 
just the energies of high-energy terms or confi gurations that ha d 
not been experimentally measured, as done bv lBadneUl (l2006ah . 

We augmented the CI expansions of the structure calcula- 
tions by including all relevant core excitations for the target ion, 
allowing for promotions out of the 45 and 4p subshells. One- 
electron additions to all target configurations were added to the 
(A^ + l)-electron ion CI expansions to account for final states for 
the Rydberg electron to radiate into the core. 

We determined the DR rate coefficients over the temperature 
range (10' - 10^)z^ K for each state within the ground configu- 
ration. Captures into Rydberg states up to « = 1000 and / - 200 
were computed, with / > 9 treated hydrogenically for all ions 
except Se^^, in which Z > 11 states were treated hydrogenically 
to allow the calculation to converge. 

3. Results 

3.1. Pliotoionization cross sections 

The photoionization and photoexcitation cross sections from 
AUTOSTRUCTURE were processed with the ADASPI and 
ADASPE codes, which respectiv ely produce the adf39 and adf38 
output files (see I Summer sl l2005l for a detailed description). To 
simplify the presentation of these results, we used the xpeppi 
post-processing cod^ to add the direct PI and photoexcitation- 
autoionization cross sections. 

PI cross sections from to 100 Rydbergs for Se"-Se^^ are 
presented as supplementary data files available at the CDS, and 
are illustrated in Fig.[r|for PI out of the ground state of each ion. 
The photoionization plus photoexcitation-autoionization cross 
sections (files labeled "se^H- _XPITOT.dat," ^ = - 5) are pro- 
vided for all states in the ground configuration of each ion. In 
those files, column 1 gives the photon energy relative to the 
ground state and column 2 the photon energy relative to the 
initial state. Columns 3 and 4 present the photo-electron en- 
ergy relative to the ground and initial states, respectively, and 
in column 5 the cross section is given in Mb. All energies are 
in Rydbergs. The data for different states are delimited by com- 
ment lines (denoted by a # sign in the first column), followed by 
three numbers indicating the energy of that level relative to the 
ground state (or the total energy, in the case of the ground state), 
the total energy of the target state, and the energy order of the 
initial state (1 for the ground state, 2 for the first excited state, 
etc.). 

We also provide direct PI cross sections without 
photoexcitation-autoionization resonances, in files labeled 

- see http://amdpp.phys. strat h.ac.uk/autos/1 
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Fig. 1. Ground state PI cross sections for Se"-Se^^, in the energy 
region near the ground state ionization threshold. 



"seq'+_XDPITOT.dat." In these files, there are two columns, the 
first of which is the photon energy relative to the ground state, 
the second the cross section in Mb. Cross sections for different 
states in the ground configuration are delimited in the same 
manner as for the se^+_XPITOT.dat files. 
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Fig. 2. Comparison of Se^^ ground state PI cross sections up to 
5 Ryd above the threshold using different CI expansions, as indi- 
cated in each of the top three panels. The bottom panel compares 
the direct PI cross sections, with photoexcitation-autoionization 
resonances removed, for the three different CI expansions. 



3.1 .1 . Estimation of uncertainties 

Special effort was made to estimate uncertainties in the PI 
cross sections, as these uncertainties affect the accuracy of 
abundances derived from spectra of astrophysical nebulae (e.g., 
ISterling et al.ll2007h . We compared the cross sections computed 
with the three CI expansions for each ion, and also varied in- 
ternal parameters or tested assumptions in our calculations (e.g., 
the values of the scaling parameters for continuum orbitals, the 
sensitivity to the gauge used, and our orthogonalization of the k- 
averaged orbitals). We also compared our results to experimental 
measurements, as described in the following subsection. 

Because we used /c-averaged relativistic wavefunctions, it 
was necessary to compute the cross sections in the velocity 
gauge even at the highest energies, where acceleration gauge 
is optimal. We compared cross sections computed in the length 
gauge at all energies to those computed solely in velocity gauge. 
We found negligible differences (<< 1%) in the cross sections 
at both low and high energies, indicating that our use of velocity 
gauge up to 100 Ryd is not a significant source of error 

For continuum orbitals, we used scaling parameters equal to 
those of the highest principal quantum number with the same 
orbital angular momentum for s, p, and d orbitals, and unity 
for higher angular momentum states. We tested this assump- 
tion by calculating PI cross sections with all continuum orbital 
scaling parameters set to unity. The direct cross sections from 
these two calculations agree to within 10% (often better) near the 



ground state threshold for all ions considered, although the res- 
onance strengths vary (especially for Se^, where the resonances 
are weaker by nearly a factor of two when the continuum scaling 
parameters are 1.0). 

For each ion, we also calculated the cross sections without 
orthogonalizing the radial orbitals. Compared to our calculations 
with orthogonalized orbitals, the direct PI cross sections agree to 
within 15% near the ground state threshold of Se**^, and to better 
than 10% for other ions. However, threshold and resonance en- 
ergies, as well as resonance strengths, are altered when orthog- 
onalization is not forced, most significantly for the neutral and 
near-neutral cases. Threshold energies tend to agree best with 
NIST values for our calculations with orthogonalized radial or- 
bitals, particularly for the neutral and single ion, in which the 
thresholds are 0.005-0.02 Ryd higher than in the non-orthogonal 
calculations. 

The most significant differences are found from comparing 
the PI cross sections calculated with different CI expansions. The 
deviations in the electronic structure for the target and (A^ -i- 1)- 
electron ions lead to disparate resonance energies and strengths, 
as can be seen in Figs. |2] and [3] for the cases of Se^^ and Se^*^. 
The sensitivity of the resonance strengths and positions to the 
CI expansion and other internal parameters illustrates the chal- 
lenge of accurately reproducing photoexcitation-autoionization 
resonances in the distorted-wave approximation. 

However, given the uncertainty in PI resonance positions 
from any type of theoretical calculation, numerical models of 
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Fig. 3. Comparison of Se^^ ground state PI cross sections up to 
5 Ryd above the threshold using different CI expansions, as indi- 
cated in each of the top three panels. The bottom panel compares 
the direct PI cross sections, with photoexcitation-autoionization 
resonances removed, for the three different CI expansions. 



ioniz ed astrophysical ne bulae often only consider direct PI 
(e.g., IVerner et al.l 1 19961) or convolve the cross section with 
Gaussians to smooth over the uncertainty in resonance positions 
(Bautista & Kallman 2001). Indeed, our goal in this work is to 
best reproduce the direct PI cross sections for low-charge Se 
ions, and to quantify their uncertainties within the distorted wave 
approximation. 

The bottom panels of Figs.|2]and[3]compare the direct ground 
state PI cross sections calculated with each of the three CI expan- 
sions. In general, the cross sections agree to within roughly 30% 
for most ions near the ground state ionization threshold. The dis- 
crepancy can be as large as 60% in the cases of Se^ and Se^^, 
but generally decreases with energy. For Se"*^, the direct cross 
sections agree to within 10% near threshold. 

Given these internal tests of our cross sections and our com- 
parisons to experimental measurements (Sec. l3.1.2t . we estimate 
our calculated direct PI cross sections to be uncertain by 30-50% 
in the energy region near the ground state ionization threshold. 



3.1.2. Comparison to experiment 

To benchmark and test the accuracy of our results, we have 
compared our (medium CI expansion) calculations to experi- 
mental absolute PI cross-section measurements performed at the 
ALS synchrotron radiation facility. Se^, Se^^, Se-'^, and Se^^ 



cross section^ were measured from the threshold of the highest- 
energy metastable state in the ground configuration up to 10 eV 
beyon d the ground state ionization threshold (Este ves et al.l 
l2010h. Experimental res ults for Se+ (iSterUng et al.l |20TT1) and 
Se-'^ jEsteves et al.ll2009 ) have been published, and we compare 
our theoretical results to experiment for those ions in this sec- 
tion (a similar comparison to the as-yet unpublished Se^^ and 
Se^^ cross sections was also performed). 

The experimen tal measurements wer e described by 
lEsteves et aU (l2009h and lSterling et al.1 (1201 ll) . and we refer the 
reader to those papers for full details. Briefly, the P I cross sec- 
tions were measured via the merged beam method (iLvon et al.l 
1986). The photoionization spectra were measured by scanning 
across photon energies, and normalized to an absolute scale 
with absolute cross-section measurements performed at discrete 
energies. Overall, uncertainties in the experimental cross 
sections were estimated to be ~30% for Se^ and ~20% for Se^^ 
(ISterling et al.ll201 itlEsteves et al.ll2009h . 

The experimental data are complicated by the fact that the 
primary ion beams were composed of an unknown admixture 
of ground and metastable states. We therefore linearly combined 
the computed cross sections of the ground configuration states to 
best reproduce the direct experimental cross sections and (when 
possible) families of resonances. 

Figs.|4]and|5]compare our calculated Se^ and Se^^ PI cross 
sections (solid lines) with experimental measurements (dotted 
lines). The computed Se^ cross section was convolved with a 
Gaussian of 28 meV FWHM to reproduce the experimental reso- 
lution. By visual comparison to the experimental data, we found 
best agreement by weighting the contribution of the ground con- 
figuration states C*S3/2, ^D3/2, ^05/2, ^Pi/2, ^P3/2 ) by (0.53, 0.15, 
0.05, 0.11, 0.16). As expected (see Sect. 13.1.11 ). our distorted- 
wave calculations do not accurately reproduce the resonance 
structure of the experimental PI cross section. In addition, the 
Fano profile near 1.7 Ryd is not reproduced by our calculations 
(R matrix techniques are necessary to reproduce resonance in- 
terference effects), and the lack of this interference causes the 
discrepancy in the direct cross section in this energy region. 
The calculated and experimental Se^ direct cross sections agree 
to within 30-50% over the range of experimental energies, ex- 
cept above 2.2 Ryd, where the experimental measurements are 
more uncertain since only a single absolute measurement was 
made above 1.98 Ryd (compared to absolute measurements at 
lower energies, where independent estimates were performed 2- 
3 times for each energy). 

The calculated Se^^ cross section was convolved with a 
20 meV FWHM Gaussian, and we found that statistical weight- 
ing of the two ground configuration states provided reasonable 
agreem ent with the experimental measurements of Esteves et al.l 
(l2009l) near threshold. Again, resonance energies and strengths 
are in poor agreement, though several general resonance fami- 
lies are seen in both the theoretical and experimental results. The 
computed direct cross section is smaller than that of experiment, 
but agrees within the noise up to 3.38 Ryd. From 3.38 to 3.9 Ryd, 
however, an apparent threshold is observed in the experimental 
data that is not reproduced in our theoretical cross section. In 
the calculated cross section, the threshold near 3.9 Ryd is due to 
photoionization into the target Se'*^ 4s4p^P term. The energy of 
this term is lower than the NIST value, and hence an inaccurate 
energy for this threshold cannot explain the experimental thresh- 



^ The Se'*^ cross section could not be measured, due to the presence 
of strong autoionizing states in the parent Se''^ ion beam that led to an 
extremely high background. 
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1.4 1.6 1.8 2.0 2.2 

Energy (Ryd) 



Fig. 4. Comparison of calculated (solid line) and experimental 
(dotted line) Se+ PI cross sections. The theoretical curve is a 
linear combination of the cross sections from the ground config- 
uration C*S3/2, ^D3/2, ^D5/2, 'P\/2, ^^3/2) Weighted by the factors 
(0.53, 0.15, 0.05, 0.11, 0.16), respectively. The theoretical cross 
section was convolved with a Gaussian of 28 meV FWHM to 
reproduce the experimental resolution. 




3.2 3.4 3.6 3.8 4.0 

Energy (Ryd) 



Fig. 5. Comparison of calculated (solid line) and experimental 
(dotted line) Se^^ PI cross sections. The theoretical curve is a 
statistically-weighted linear combination of the cross sections 
from the two ground configuration states. The theoretical cross 
section was convolved with a Gaussian of 20 meV FWHM to 
reproduce the experimental resolution. 



old at 3.38 Ryd. Given that 4s 4p is the lowest-energy excited 
term of Se"*"^ and that it is unlikely the primary Se^^ ion beam 
contained metastable states with excitation energies of more than 
one Rydberg, this experimental feature is puzzling. However, re- 
cent revisions to the analysis of the experimental Se^^ data indi- 
cate that the apparent threshold at 3.38 Ryd is due to a cluster of 
poorly-resolved resonances, as well as an error in the normaliza- 
tion of the data to an absolute scale at energies above 3.38 Ryd 
(D. A. Esteves, private communication). 

Comparisons to the as-yet unpublished Se^^ and Se^^ cross 
sections (not shown in this paper) are more encouraging, at least 



in terms of the direct cross sections. In both cases, our calcula- 
tions agree within the noise of the experimental data. 

3.2. Radiative recombination rate coefficients 

RR rate coefficients were calculated from the direct photoion- 
ization cross sections using the ADASRR code. The resulting 
adf48 output files (availa ble at the CDS) are similar in stiTicture 
to adp9 files for DR (see[S ummer 120051 for a detailed descrip- 
tion). In these files, partial RR rate coefficients to states with 
« < 8 are fully /-resolved, and are bundled over the final outer 
quantum numbers as bundled-«/ for « < 10, and bundled-n for 
n < 1000 (see .BadnelL2006b ). At the end of each adf48 file, to- 
tal RR rate coefficients are written as a function of temperature 
for each target ground configuration state. 

We fit the rate coefficients with the analytical form 
(IVemer & Ferlandlll996tlB"adnellll2006bh 

ffRR(r) = A X [ ^/f/To{l + ^lf|n,)'-''Hl + ^lfWl)'^''«Y ' 
where 

Bo = B + Ce\p{-T2/T). (2) 

In this equation, To \ 2 are in temperature units, A and the rate co- 
efficient aRR(r) have units of cm^ s and B and C are dimen- 
sionless parameters. For several low-charge states, the RR rate 
coefficients exhibit large "bumps" at high temperature. These 
features are caused by electron capture into low principal quan- 
tum number states, in which the nucleus is not as effectively 
screened as for higher states. Equations [1] and |2] are not able to 
reproduce the size of the bumps, and become inaccurate at high 
temperatures. For these Se^-Se^^ target states, we therefore fit 
the high-temperature RR rate coefficients with the function 

1 " 

aRR(T) = ^ c,exp(-£',/r), (3) 

/-I 

where T and E, are in temperature units (K), the rate coefficient 
is in cm-' s"', c,- in cm-' s"' K^^^, and n ranges from 5 to 7 de- 
pending on the ion. This ana lytical function is typic ally used to 
fit DR rate coefficients (e.g., IZatasrinnv et al.l[2003h . as done in 
Sect. [33] 

The fit coefficients were determined using a non-linear least- 
squares fit algorithm, generally leading to fit accuracies within 
5% for the low-temperature fits, and <2% at high temperatures. 
However, the low-temperature fits are only accurate to within 
7.5% for Se^ (level 1) and Se-'^ (level 1), and the accuracies 
of the high-temperature fits for Se^ are 3% (levels 2-4) to 4% 
(level 5). Moreover, Equations [1] and |2] are unable to reproduce 
the rate coefficient for RR onto level 5 of Se^ at the lowest tem- 
peratures (< 20 K), where they deviate from the calculated rate 
coefficients by up to 20%. Aside from this exception, the fits 
have the correct asymptotic forms outside of the temperature 
range (10' - 10^ )z^ K (specifically. Equations [T] and [2] describe 
the asymptotic behavior as T — > K, and Equation [3] correctly 
describes the rate coefficient at the high-temperature limit). 

Table |5] displays the fit coefficients from equations [T] and 
|2] as well as the maximum temperature T^ax for which the fits 
are valid to within the stated accuracies. If Tmax is not given in 
Table|5] then equations[T]and|2]are accurate over the full temper- 
ature range (10' - 10^)z^ K and coiTectly describe the asymptotic 
behavior of the rate coefficients at the high-temperature limit. 
Table |6]provides fit coefficients for the RR rate coefficients from 
T^max to lO^z^ K. In Fig.[6l we display the ground state RR rate 
coefficients for low-charge Se ions as a function of temperature. 
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Table 5. Fit coefficients for radiative recombination rate coefficients at low temperature (see Equations [1] and |2). r^ax represents 
the maximum temperature for which these fits are valid to within the stated accuracies (5%, unless noted). For recombination at 
temperatures exceeding Tmax, the fits from Equation[3]and coefficients listed in Table |6] should be used. If no Tmax is given, then the 
fits are accurate over the entire temperature range (10' - 10^)z^ K. The notation x(y) denotes x x W. 



Target 




A 




B 


To 


Ti 


C 


r. 


T 

max 


Ion" 


Level 


(cm^ s 


-') 




(K) 


(K) 




(K) 


(K) 


Se+ 


1'^ 


7.395(- 


■10) 


0.7288 


7.273(-2) 


3.495(6) 


0.3306 


2.160(6) 






2 


2.808(- 


11) 


0.4980 


3.761 


3.415(129) 


0.3365 


4.621(1) 


2.0(3) 




3 


2.747(- 


■11) 


0.4818 


3.926 


2.747(111) 


0.3515 


4.451(1) 


2.0(3) 




4'- 


1.439(- 


11) 


0.0535 


2.036(1) 


6.815(130) 


0.7152 


4.065(1) 


2.0(3) 




5" 


3.719(- 


■11) 


0.4765 


3.786 


3.019(101) 


0.2935 


6.541(1) 


1.0(3) 


Se2+ 


1 


6.253(- 


11) 


0.6744 


6.081(1) 


5.440(5) 


0.8633 


2.232(6) 


8.0(4) 




2 


3.804(- 


12) 


0.5422 


7.174(3) 


3.930(6) 


1.3093 


1.915(5) 


8.0(4) 




3 


2.977(- 


■11) 


1.0695 


3.398(1) 


1.907(5) 


0.9534 


2.179(6) 


8.0(4) 




4 


1.358(- 


11) 


1.1128 


7.719(1) 


2.460(5) 


1.0094 


2.152(6) 


8.0(4) 




5 


2.109(- 


11) 


1.0811 


2.623(1) 


2.879(5) 


0.8678 


2.142(6) 


8.0(4) 


Se3+ 


1<- 


1.496(- 


11) 


0.0000 


3.270(3) 


9.004(7) 


0.7282 


1.394(6) 






2 


5.353(- 


11) 


0.9825 


5.028(1) 


6.721(5) 


1.0127 


4.915(5) 


9.0(5) 


Se-^^ 


1 


2.664(- 


■11) 


0.0000 


3.018(3) 


1.142(8) 


0.6892 


1.136(6) 




Se5+ 


1 


5.822(- 


11) 


0.1471 


1.562(3) 


1.905(8) 


0.5119 


1.031(6) 




Se*+ 


1 


1.293(- 


■10) 


0.3096 


8.107(2) 


2.554(8) 


0.3114 


1.278(6) 





Notes. Note that the ion and level numbers correspond to the target ion (i.e., before recombination). The fit is accurate to within 7.5%. The 
fit is accurate to within 5% except at the lowest temperature (10 K), where it is accurate to within 7%. '''' The fit is accurate to within 5% except at 
the lowest temperatures (below 20 K), where the fit agreement can be as poor as 20%. For this reason, the fit to the RR rate coefficent for this state 
should not be considered valid below 20 K. '""^ The fit is accurate to within 7%. 

Table 6. Fit coefficients for radiative recombination rate coefficients at high temperatures, valid from r^ax (see TablejSj to lO^z^ K. 
See Equation [3] Coefficients Ci are in cm-^ s"' K^^^, and Ei are in K. States whose RR rate coefficients are accurately described by 
Equations [T] and 12] (with coefficients given in Table|5]l over the full temperature range (10' - 10^)z^ K are not listed here. The fits 
are accurate to within 3% for Se^ levels 2^, 4% for Se^ level 5, and <2% otherwise. The notation x(y) denotes x x 10-'. 



Target 



Ion" 


Level 




'2 


C3 


Q 


cs 


C6 




Se+ 


2 


6.380(-8) 


3.641(-7) 


1.412(-6) 


2.856(-6) 


6.461(-6) 


1.067(-4) 


4.751(-4) 




3 


6.425(-8) 


3.638(-7) 


1.403(-6) 


2.81 1(-6) 


6.138(-6) 


1.018(-4) 


4.729(-4) 




4 


7.502(-8) 


4.203(-7) 


1.429(-6) 


2.582(-6) 


5.538(-6) 


8.929(-5) 


3.948(-4) 




5 


3.371(-8) 


1.857(-7) 


8.251(-7) 


2.939(-6) 


5.302(-6) 


8.719(-5) 


3.936(-4) 




1 


2.061(-6) 


8.034(-6) 


1.447(-5) 


2.769(-4) 


1.237(-3) 


2.609(-3) 






2 


2.232(-6) 


8.55 1(-6) 


1.552(-5) 


2.780(-4) 


1.240(-3) 


2.586(-3) 






3 


1.018(-5) 


1.613(-5) 


2.749(-4) 


1.236(-3) 


2.612(-3) 








4 


1.195(-5) 


7.588(-5) 


6.316(-4) 


1.965(-3) 


4.432(-3) 








5 


5.306(-6) 


7.409(-6) 


7.781(-5) 


5.502(-4) 


1.469(-3) 


2.035(-3) 




Se'-^ 


2 


4.289(-5) 


1.638(-3) 


1.161(-2) 


5.245(-3) 


2.772(-2) 










El 


E2 


E3 


E4 


Es 


Ee 


El 


Se+ 


2 


6.920(2) 


6.054(3) 


2.936(4) 


1.100(5) 


6.510(5) 


3.218(6) 


1.207(7) 




3 


6.990(2) 


6.080(3) 


2.926(4) 


1.087(5) 


6.306(5) 


3.137(6) 


1.190(7) 




4 


6.984(2) 


6.165(3) 


2.893(4) 


1.029(5) 


6.244(5) 


3.151(6) 


1.194(7) 




5 


3.045(2) 


2.838(3) 


1.405(4) 


6.763(4) 


5.084(5) 


3.075(6) 


1.182(7) 




1 


3.108(3) 


4.135(4) 


3.780(5) 


2.581(6) 


1.008(7) 


3.998(7) 






2 


4.231(3) 


4.997(4) 


4.748(5) 


2.637(6) 


1.015(7) 


4.037(7) 






3 


3.628(4) 


4.849(5) 


2.644(6) 


1.019(7) 


4.124(7) 








4 


6.274(4) 


1.334(6) 


5.185(6) 


2.057(7) 


1.135(8) 








5 


2.145(4) 


1.718(5) 


1.457(6) 


5.086(6) 


1.786(7) 


6.142(7) 




Se3+ 


2 


6.867(4) 


1.119(6) 


3.421(6) 


2.101(7) 


9.183(7) 







Notes. Note that the ion and level numbers correspond to the target ion (i.e., before recombination). 



3.2.1. Estimation of uncertainties 



We tested the sensitivity of the calculated RR rate coefficients 
to the use of different CI expansions and other internal parame- 
ters in our calculations. The uncertainties that we discuss were 
estimated at temperatures near 10"* K, similar to electron temper- 
atures in photoionized nebulae. 



As is the case for PI, the RR rate coefficients are most sen- 
sitive to the CI expansion used. The rate coefficients calculated 
with the small and large configuration sets (Table [1} differ from 
the presented rate coefficients (from the medium CI expansion) 
typically by less than 10% near lO** K. The sole exception is that 
of Se^, for which the deviation from the presented rate coeffi- 
cients is 25-30%. 
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10^ 10^ 10=^ 10^ 10^ 10'^ 10'' lO*' 10^ 
Temperature (K) 

Fig. 6. Radiative recombination rate coefficients as a function of 
temperature for ground state target ions Se+ (solid curve), Se^^ 
(dotted curve), Se^*^ (short-dashed cui-ve), Se'*^ (dash-dot-dash 
curve), Se^^ (long-dashed cui-ve), and Se^^ (dash-dot-dot-dot- 
dash curve). 

Other tested parameters produced smaller changes in the RR 
rate coefficients. For example, to test the sensitivity of forcing 
velocity gauge at high energies, we computed the rate coeffi- 
cients from PI cross sections cut off at the maximum velocity 
gauge energy (beyond which velocity gauge must be forced in 
AUTOSTRUCTURE). The RR rate coefficients were unaffected 
except at the highest temperatures (typically a few times 10^- 
10^ K). This indicates that uncertainties associated to our use 
of velocity gauge at high energies are negligible for the com- 
puted RR rate coefficients at photoionized plasma temperatures, 
and become important only above 10^ K for singly- and doubly- 
charged Se, and above 10^ K for higher charge states. 

The calculations were performed using 3 interpolation ener- 
gies per decade. We tested the effects of using a finer interpo- 
lation mesh of five energies per decade, and found differences 
of less than 1 % in the RR rate coefficients for all of the investi- 
gated Se ions. Similarly, the use of target ion scaling parameters 
versus those of the + 1 -electron ion produced discrepancies in 
the rate coefficients larger than 1% only in the case of Se^ and 
Se^^, where the differences are 3-7% (depending on the state) 
near 10"* K. Finally, we computed RR rate coefficients without 
Schmidt orthogonalizing the radial orbitals, and found that the 
rate coefficient differed by a maximum of 1-3% (for Se^ and 
Se^^) compared to our orthogonalized calculations. 

We conclude that our computed RR rate coefficients have 
internal uncertainties of less than 10% with the exception of Se^, 
where the uncertainty is 30-40%. 

3.3. Dielectronic recombination rate coefficients 

DR rate coefficients were computed from the resonant photoion- 
ization cross sections with the ADASD R cod e, and the adf09 
output files are available at the CDS (see [Summ ers 2005, for a 
detailed description of the structure of these files). The rate coef- 
ficients were fitted with functions of the form Equation |3] with n 
ranging between 5 and 7 (see Table|7]i. The fits for most ions are 
accurate to within 5 % over the temperature range ( 1 ' - 1 0^ )z^ K, 




Temperature (K) 

Fig. 7. Dielectronic recombination rate coefficients as a function 
of temperature for ground state target ions Se^ (solid curve), 
Se^+ (dotted curve), Se^+ (short-dashed curve), Se^+ (dash-dot- 
dash curve), Se^^ (long-dashed curve), and Se^+ (dash-dot-dot- 
dot-dash curve). 

and reproduce the correct asymptotic behavior outside of this 
temperature range. The fitting algorithm was not able to accu- 
rately reproduce the rate coefficient behavior over the full tem- 
perature range for DR onto level 4 (the third excited state) of the 
Se^^ target. That fit is not valid below 80 K, and is only accu- 
rate to within 11% near 2 000 and 40000 K (but to better than 
2% at other temperatures). DR onto level 2 of Se^ is also fit less 
accurately, to within 6%. 

Fig. |7] illustrates the DR rate coefficients over the temper- 
ature range (10' - 10^)z^ K. This figure reveals that while the 
rate coefficients exhibit similar behavior at high temperatures 
(> 10^ K), there are marked differences in their magnitudes and 
behavior at lower temperatures due to the differing resonance 
structure near their ionization thresholds. 

Fig. [8] compares the RR, DR, and total recombination rate 
coefficients for each of the first six Se ions. From this plot, it 
is apparent that DR dominates RR for each ion at lO'* K, with 
rate coefficients larger by as much as two orders of magnitude 
(and occasionally even more at other temperatures). The RR rate 
coefficient becomes comparable to that of DR at temperatures 
higher than 10^ K for some ions, though this is likely due to the 
fact that we did not consider DR from An > core excitations, 
which become important at high temperatures. RR dominates 
DR only at low temperatures (less than 1 000-5 000 K) for Se+, 
Se^^, and Se*^. This highlights the importance of obtaining ac- 
curate DR rate coefficients, as it is the principal recombination 
process for low-charge Se ions (and indeed for many other heavy 
elements, as noted for ex ample by |Badnell 20 061 lAltun et al.l 
120071 : iNikolic et al]|2010l) in photoionized nebulae. 

3.3.1. Estimation of uncertainties 

At photoionized plasma temperatures, DR proceeds via elec- 
tron capture into low-lying autoionizing states. Unfortunately, 
the energies of these states have not been determined for 
the vast majority of atomic ions beyond the second row of 
the Periodic Table. Inaccuracies in the calculated energies of 
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Table 7. Fit coefficients for dielectronic recombination rate coefficients; see Equation[3] Coefficients q are in cm^ s"' K^^^, and E; 
are in K. The fits are accurate to within 5% over the temperature range (10^ - 10^)z^ K, unless otherwise noted. The notation x(y) 
denotes x x 10' . 



Target 



Ion" 


Level 




Cl 


C3 


C4 


C5 




c? 


Se+ 


1 


5.895(-8) 


1.516(-7) 


2.154(-6) 


4.509(-6) 


7.933(-4) 


2.944(-3) 


2.944(-3) 




2* 


4.849(-9) 


4.750(-8) 


1.915(-7) 


1.839(-6) 


1.633(-4) 


2.211 (-3) 


7.745(-6) 




3 


2.580(-8) 


4.667(-8) 


2.042(-6) 


2.148(-3) 


5.616(-4) 








4 


2.086(-8) 


2.435(-7) 


2.120(-6) 


3.332(-7) 


3.652(-6) 


1.713(-3) 


1.929(-4) 




5 


2.547(-7) 


1.297(-7) 


4. 131 (-4) 


2.075(-3) 


2.433(-6) 






Se2+ 


1 


4.841(-6) 


4.993(-5) 


3.566(-5) 


9.558(-5) 


9.458(-3) 


1.988(-5) 






2 


1.830(-6) 


5.312(-6) 


1.432(-5) 


2.280(-5) 


5.604(-3) 


2.176(-3) 






3 


1.686(-8) 


2.146(-7) 


5.827(-7) 


2.436(-6) 


1.391(-4) 


6.526(-3) 






4' 


3.875(-8) 


1.641(-6) 


3.522(-3) 


3.732(-3) 


-3.333(-4) 








5 


7.246(-13) 


9.784(-9) 


3.270(-6) 


2.470(-5) 


2.923(-4) 


1.598(-3) 


1.609(-3) 


Se3+ 


1 


1.268(-6) 


2.651(-6) 


1.386(-5) 


1.166(-4) 


3.731(-3) 


6.768(-3) 


6.499(-3) 




2 


2.286(-6) 


3.420(-6) 


4.948(-5) 


4.206(-4) 


1.030(-2) 






Se^+ 


1 


1.505(-5) 


3.440(-5) 


2.571(-4) 


2.806(-3) 


2.952(-2) 






Se^-*- 


1 


1.107(-5) 


2.896(-5) 


8.937(-5) 


6.877(-5) 


1.457(-2) 


1.213(-2) 




Se«+ 


1 


2.771(-6) 


3.733(-5) 


1.459(-4) 


7.999(-4) 


5.106(-3) 


1.052(-2) 












Ei 


£4 


Es 


Ee 


El 


Se+ 


1 


1.036(3) 


1.963(3) 


4.529(3) 


1.464(4) 


1.217(5) 


1.754(5) 


1.772(5) 




2" 


8.714(1) 


2.241(2) 


7.538(2) 


1.046(4) 


9.325(4) 


1.525(5) 


1.764(6) 




3 


9.716(1) 


1.663(3) 


1.077(4) 


1.330(5) 


2.180(5) 








4 


2.655(2) 


3.997(2) 


9.086(2) 


1.797(3) 


3.584(4) 


1.336(5) 


2.199(5) 




5 


8.821(2) 


2.098(3) 


9.386(4) 


1.423(5) 


8.491(5) 






Se2+ 


1 


1.289(2) 


5.862(2) 


1.935(3) 


7.117(4) 


1.957(5) 


1.719(6) 






2 


6.188(1) 


3.923(2) 


2.220(3) 


3.978(4) 


1.797(5) 


2.263(5) 






3 


4.308(2) 


9.211(2) 


3.360(3) 


1.096(4) 


8.181(4) 


1.899(5) 






4" 


4.521(3) 


1.678(4) 


1.429(5) 


2.330(5) 


5.738(5) 








5 


1.158(3) 


5.031(3) 


1.522(4) 


3.687(4) 


9.171(4) 


1.416(5) 


1.657(5) 


Se3+ 


1 


7.565(2) 


1.959(3) 


6.654(3) 


3.054(4) 


1.518(5) 


2.134(5) 


2.156(5) 




2 


6.470(2) 


3.122(3) 


1.519(4) 


6.767(4) 


1.908(5) 






Se^+ 


1 


5.954(2) 


3.483(3) 


2.477(4) 


1.053(5) 


1.886(5) 






Se5+ 


1 


3.417(2) 


2.862(3) 


1.556(4) 


4.233(4) 


1.420(5) 


1.941(5) 




Se**+ 


1 


2.567(3) 


8.368(3) 


2.392(4) 


1.260(5) 


5.466(5) 


1.186(6) 





Notes. Note that the ion and level numbers correspond to the target ion (i.e., before recombination). The fit is accurate to within 6%. The 
rate coefficient for DR onto level 4 of the Se^+ target could not be fit to accuracies as high as other levels and ions, and is not valid below 80 K. 
The fit is only accurate to within 11% near 2 000 and 40 000 K, but to within 2% at all other temperatures. 



these near-threshold autoionizing resonances are thus the chief 
source of uncertainty in computed low-temperature DR rate co- 
efficients (see, 



Savin et all \l99^ iFerlandl I200I . While 



IRobicheaux et al.N2oiol) showed that thermal broadening of the 
continuum threshold in finite density plasmas alleviates to some 
extent the problem of identifying which resonances lie just above 
the threshold (and can contribute to DR) and which do not, this 
effect is not likely to be sufficiently large to mitigate the compar- 
atively large uncertainties in resonance positions resulting from 
theoretical calculations. 

To estimate the uncertainties in our calculated DR rate coef- 
ficients due to the unknown energies of near-threshold autoion- 
izing resonances, we shifted the continuum threshold by the 
amount of the largest discrepancy in calculated versus experi- 
mental target energies, with the direction of the shift determined 
by whether the computed energies were too large or too small 
relative to experiment. For DR onto the ground state of target 
ions, this resulted in discrepancies in the rate coefficient near 
lO-* K of 15-30% for Se^+ and Se'*+, a factor of two for Se^^, a 
factor of 4-5 for Se^^, and 1-2 orders of magnitude for Se^^ and 
Se^ (respectively). 

We also compared the DR rate coefficients computed with 
the small and medium configuration sets (due to the compu- 
tationally intensive nature of the DR calculations, we did not 



attempt to compute rate coefficients for the large configuration 
sets). These uncertainties were much smaller than obtained from 
shifting the continuum threshold, leading to differences of 5- 
30% near 10"^K with the exception of the Se^ target (factor of 4 
difference). 

Clearly, the uncertainties in the DRrate coefficients are much 
larger than those for PI or RR. Because DR dominates RR at 
photoionized plasma temperatures for these Se ions (and in- 
deed many other heavy atomic ions), this remains an important 
problem for modeling photoionized astrophysical nebulae. Since 
the energies of near-threshold autoionizing resonances have not 
been experimentally determined and are extremely difficult to 
accurately calculate even with the most sophisticated theoretical 
methods, the solution to this important problem will be challeng- 
ing to achieve. 

Perhaps the most reliable method of obtaining highly ac- 
curate DR rate coefficients is electron-ion beam mergin g tech- 
nique s with heavy-ion storage rings s uch as the TSR (Savin et al.l 
1 19991 l2006t ISchippers et 5] l201 ih . Unfortunately, it is ex- 
tremely challenging to use storage ring experiments to measure 
DR rate coefficients of low charge-to-mass species such as the 
Se ions investigated in this paper. The lowest charge-to-mass ion 
whose recombination rate coefficient h as been experimentall y 
determined is Sc^+, with q/m - 1/15 dSchippers et all l2002l) . 
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Fig. 8. Comparison of the radiative (dotted curves), dielectronic 
(dashed curves), and total recombination rate coefficients (solid 
curves) for each of the first six Se ions. 



The maximum bending power of the dipole magnets on cur- 
rent storage ring facilities is the Umiting factor for low qjm 
species, which cannot be stored at sufficiently high energies. If 
the ion energy is too low, then electron capture during collisions 
with residual gas molecules produces a very high background 
that prevents rate coe fficient measurements (S. Schippers, pri- 
vate communication; ISchippers et al.ll2b02h . Thus, experimen- 
tal DR measurements of low-charge Se ions are extremely dif- 
ficult if not impossible with existing storage ring facilities. 
However, the forthcoming Cryogenic Storage Ring (CSR) fa- 
cility in Heidelberg will use electrostatic deflectors in place of 
bending magnets to store the ions, and the residual gas pres - 
sure will be much smaller than in the TSR (Wolf et aDl2006h . 
The CSR will enable low-temperature DR rate coefficients to be 
measured for many low qjm species, including those of trans- 
iron elements. 



4. Concluding remarks 

We have presented MCBP calculations of the electronic struc- 
ture, distoited-wave PI cross sections (including both direct and 
resonant contributions), and RR and DR rate coefficients for the 
first six ions of the trans-iron element Se. All of the PI, RR, and 
DR data are available at the CDS, and analytical fits to the RR 
and DR rate coefficients are presented in Tables |5]|6] and|7] 

These atomic data are critical for efforts to determine 
the abundance of Se in photoionized astrophysical nebu- 



lae. Se has been de t ected in dozens of planetary nebulae 
(ISterling & DinersteinI 12008! ) and other astrophysical nebulae 
(see Sect. [1) , and is a potentially important tracer of «-capture 
nucleosynthesis and the chemical evolution of trans-iron ele- 
ments. 

In our calculations, we tested the sensitivity of our results 
to the use of different configuration expansions and to other in- 
ternal parameters of the AUTOSTRUCTURE code, in order to 
provide estimates of the (internal) uncertainties of the presented 
atomic data. In the case of PI, we also com pared our cross sec- 
tions to recent exper imental measurements (lEsteves et al.ll2009l 
2010; Sterling et al. 2011) to gauge the accuracy of our results. 
We find that the direct PI cross sections for these Se ions are un- 
certain by 30-50%, while the RR rate coefficients are uncertain 
by <10% except for Se+ (-30-40%). Our DR rate coefficients 
exhibit larger uncertainties (from 15-30% up to 1-2 orders of 
magnitude) due to the unknown energies of low-lying autoion- 
izing states. The importance of uncertainties in DR rate coef- 
ficients is highlighted by the dominance of DR over RR near 
10** K, the typical temperature of photoionized nebulae such as 
PNe. 

These uncertainties are significant, particularly for the near- 
neutral cases. Low-charge Se ions are complex systems for 
which no comprehensive photoionization or recombination data 
existed prior to our study. Thus one of our goals in addition 
to producing these data was to provide realistic internal uncer- 
tainties, whose effect on nebular abundance determinations can 
be quantified through numerical simulations of ionized nebulae. 
Such an investigation will help to target the ionic systems and 
atomic processes that require further analysis. 

This paper is the first in a series to present atomic data for 
the photoionization and recombination properties of trans -iron 
elements. In future papers, we will present similar data for Kr 
and Xe, as well as charge exchange rate coefficients for low- 
charge ions of several «-capture elements. The ultimate goal of 
this study is to produce atomic data for n-capture elements that 
are suitable to be incorporated into photoionization codes that 
numerically simulate the ionization and thermal structure of as- 
trophysical nebulae. Photoionization codes can be used to derive 
analytical ionization corrections needed to accurately estimate 
the abundances of unobserved Se, Kr, and Xe ions in ionized 
nebulae - and hence to determine their elemental abundances. 
Our efforts will thus enable much more accurate and robust 
abundance determinations of trans-iron elements in astrophys- 
ical nebulae than was previously possible. 
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Table 1. Configuration-interaction expansions used for each Se ion. 

Ion Config. Set" Configurations 
~~S^ Small (10) 4sMp\ 4s4p\ 4p^ 

Medium (326) 4s^ 4p\ 4s^ 4p^ 4d, 4s^ 4p^ 4f, 48^ 4p^ 5s, 4s- 4p'' 5p, 

4sMpMd5p, 4s4p^ 4s4p'*5p, 4s4p'5p^ 4p* 
Large (785) 4sMp\ 4sMpMd, 4s2 4pMf, 4s2 4p3 5s, 4sMp3 5p, 

4s2 4p3 5d, 4s2 4p2 4d2, 4s2 4p2 4d5p, 4s4p^ 4s4p''4d, 
4s4p''5s, 4s4p**5p, 4s4pMd5s, 4s4p3 5p^ 4p^ 

3d^ 4s^ 4p^ 

~~Se^ Small (15) 4?4p7^s4p'7^p 

Medium (189) 4sMp^ 4sMpMd, 4s2 4p2 5s, 4s2 4p2 5p, 4s-4p4d^ 

4s4p\ 4s4p3 4d, 4p^ 
Large (886) 4sMp^ 4sMpMd, 4sMp2 5s, 4s2 4p2 5p, 4sMp2 5d, 
4sMp4d2, 4s2 4p4d5s, 4sMp4d5d, 4s4p\ 4s4pMd, 
4s4pMd^ 4p^ 4pMd, 4p-*5s 



Se-+ Small (32) 


4sMp^ 4sMp4d, 4s4p^ 4p'* 


Medium (135) 


4s^4p2, 4s2 4p4d, 4s 4p^ 4s4p^4d, 4p**, 




4pMd2, 3d''4sMp3 


Large (765) 


4sMp^ 4sMp4d, 4sMp5s, 4s^4p5p, 4s^4p5d, 




4sMd2, 4sMd5d, 4s4p3, 4s4p2 4d, 4s4p2 4f, 




4s4p4d2, 4p^ 4pMd, 4p3 5s, 3d''4sMp2 4d, 




3d' 4sMp4d5s 



Se^+ Small (40) 


4s^4p, 4s^4d, 4s 4p'', 4s4p4d, 4p^ 


Medium (252) 


4s^4p, 4s^4d, 4s 4p^, 4s4p4d, 4s 4d^, 




4p3, 4pMd, 3d''4s2 4p4d, 3d''4s4p3 


Large (546) 


4s2 4p, 4s2 4d, 4s2 5s, 4s 4p^, 4s4p4d, 




4s4p4f, 4s4p5s, 4s4p5p, 4s4p5d, 4s 4d^, 




4s4d5s, 4s4d5p, 4p\ 4p2 4d, 4p- 5p, 




4p4d2, 4p4d5s, 3d''4sMp4d, 3d''4sMd2, 3d'4sMd5s, 




3d''4s4p3 



Se*+ Small (10) 4s2, 4s4p, 4p2 

Medium (47) 4s^ 4s 4p, 4s 4d, 4s 5s, 4p-, 

4p4d, 4d5d, 58^ 
Large (891) 4s^ 4s4p, 4s4d, 4s5s, 4s5p, 

4s 5d, 4p^, 4p4d, 4p5s, 4p5p, 

4d^ 4d5s, 4d5d, 5s^ 5p^ 

SdMsMd, 3d''4s2 5d, 3d''4s4p2, 3d''4s4p4d, 3d''4s4p4f, 
3d'4s4f5p, 3d'' 4s 5s 5d 



Se"+ Small (3) 


4s, 4p 


Medium (5) 


4s, 4p, 4d 


Large (745) 


4s, 4p, 4d, 5s, 5p, 




3d''4s4d, 3d''4s5d, 3d''4p4d, 3d''4p4f, 3d''4p5d, 




3d''4d-, 3d''4d5s, 3d''4d5p, 3d''4d5d, 3d''5s5p 


Se'-^ Small (35)" 


3d'", 3d'' 4s, 3d'Mp, 3d'' 4d 


Medium (35)* 


3d'", 3d'' 4s, 3dMp, 3d'' 4d 


Large (57) 


3d"', 3d'' 4s, 3d''4p, 3d'' 4d,3d' 5s, 




3d''5d 



Notes. After the name of each configuration set, the number of levels is listed in parentheses. Although the same CI expansion was utilized 
for the small and medium sets of Se^^, the scaling parameters of the medium set were optimized differently: the configurations 3d' 5d, 3d' 6d, and 
3d' 7d were included in the configuration list to determine scaling parameters, but removed in the final structure calculation. 
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Table 2. Radial scaling parameters used for CI expansions of each ion 
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Table 2. Continued. 



Orbital 
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Medium 


Large 
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Table 3. Comparison of selected calculated and experimental energies (in Rydbergs) 
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0.5462 


0.5764 


0.5949 


21 








0.5478 


0.5776 


0.5954 


22 


45^ 4p3 5s 


^D2 




0.5569 


0.5907 


0.5634 


23 




^D, 




0.5571 


0.5902 


0.5621 


24 








0.5584 


0.5925 


0.5672 


25 




'D2 




0.5701 


0.6134 


0.5785 


I.P. 






0.6377 


0.6558 


0.6590 


0.7168 




Index 


Config. 


Term 


Small 


Medium 


Large 


NIST 


1 


4s- 4p^' 


*S3/2 
'D3/2 
'D5/2 


0.0000 


0.0000 


0.0000 


0.0000 


2 




0.1527 


0.1470 


0.1413 


0.1200 


3 




0.1581 


0.1520 


0.1457 


0.1256 


4 




'Pl/2 


0.2204 


0.2218 


0.2104 


0.2099 


5 




^P3/2 


0.2267 


0.2277 


0.2156 


0.2177 


6 


4s 4p** 


**P5/2 


0.7851 


0.7569 


0.7460 


0.7643 


7 




**P3/2 


0.8001 


0.7708 


0.7582 


0.7799 


8 




**Pl/2 


0.8081 


0.7783 


0.7647 


0.7877 


9 


4s2 4p2 5s 


**Pl/2 




0.9045 


0.8912 


0.8682 


10 




**P3/2 




0.9145 


0.9002 


0.8817 


11 




''P5/2 




0.9289 


0.9129 


0.8992 


12 




-Pl/2 




0.9433 


0.9254 


0.9012 


13 




-P3/2 




0.9523 


0.9276 


0.9236 


14 




'Dj/2 




1.0396 


1.0192 


0.9874 


15 




2D3/2 




1.0399 


1.0187 


0.9883 


16 


4s^ 4p- 5p 


"Si/2 




1.0541 


1.0644 


1.0699 


IP. 






1.5128 


1.5242 


1.5388 


1.5574 


Se^+ 


Index 


Config. 


Term 


Small 


Medium 


Large 


NIST 


1 


4s- 4p^ 


^Po 


0.0000 


0.0000 


0.0000 


0.0000 


2 




3P, 


0.0135 


0.0142 


0.0158 


0.0000 


3 




3P2 


0.0315 


0.0334 


0.0362 


0.0000 


4 




'D2 


0.1383 


0.1406 


0.1407 


0.1188 


5 




'So 


0.2424 


0.2710 


0.2595 


0.2591 


6 


4s 4p^ 


^D, 


0.7935 


0.8376 


0.8308 


0.8301 


7 




^D2 


0.7937 


0.8378 


0.8313 


0.8450 


8 




'D3 


0.7967 


0.8412 


0.8353 


0.8798 


9 




-^P2 


0.9267 


0.9727 


0.9654 


0.9706 


10 




^^Po 


0.9274 


0.9735 


0.9675 


0.9703 


11 




-^Pl 


0.9275 


0.9737 


0.9676 


0.9713 


12 


4s2 4p4d 


'02" 


1.0250 


1.0685 


1.0262 


1.2685 


13 




^^F2 


1.1291 


1.1691 


1.1314 


1.1304 


14 




3F3 


1.1390 


1.1795 


1.1424 


1.1419 


15 




'F4 


1.1543 


1.1962 


1.1602 


1.1610 


16 




'Pi 


1.2586 


1.3019 


1.3215 


1.2479 


17 




^P2 


1.2984 


1.3365 


1.3042 


1.3004 


18 




3P, 


1.3067 


1.3453 


1.3310 


1.3009 


19 




3Po 


1.3156 


1.3559 


1.3283 


1.2941 


20 






1.3292 


1.3673 


1.3054 


1.2816 


21 




^^D3 


1.3331 


1.3701 


1.3255 


1.2969 


22 




^D2 


1.3338 


1.3718 


1.3328 


1.2704 



N, 
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Table 3. Continued. 



23 




'F3 


1.4080 


1.4427 


1.3855 


1.3548 


24 


4s4p3 




1.4366 


1.4801 


1.4289 


1.0258 


IP. 






2.2322 


2.2591 


2.2585 


2.2653 


Se^+ 


Index 


Config. 


Term 


Small 


Medium 


Large 


NIST 


1 


4s'^4p 


'Pl/2 


0.0000 


0.0000 


0.0000 


0.0000 


2 




^P3/2 


0.0359 


0.0408 


0.0399 


0.0399 


3 


4s4p2 




0.6690 


0.6954 


0.7242 


0.7235 


4 




''P3/2 


0.6826 


0.7110 


0.7394 


0.7379 


5 




''P5/2 


0.7025 


0.7332 


0.7612 


0.7617 


6 




2D3/2 


0.9631 


0.9645 


0.9848 


0.9496 


7 




^05/2 


0.9649 


0.9675 


0.9878 


0.9542 


8 






1.2240 


1.2380 


1.2570 


1.1736 


9 




^Pl/2 


1.3662 


1.3312 


1.3517 


1.2405 


10 




^P3/2 


1.3883 


1.3537 


1.3735 


1.2608 


11 


4s2 4d 


^D3/2 


1.4846 


1.4520 


1.4404 


1.3962 


12 




'D5/2 


1.4868 


1.4542 


1.4425 


1.3998 


13 


4p3 


^D3/2 


1.7948 


1.8299 


1.8413 


1.7288 


14 




-D5/2 


1.7999 


1.8360 


1.8473 


1.7948 


16 




'P3/2 


2.0512 


2.0962 


2.0974 


1.9987 


17 




-Pl/2 


2.0536 


2.0988 


2.0977 


2.0644 


LP. 




3.0904 


3.1181 


3.1140 


3.1564 


86*+ 


Index 


Config. 


Term 


Small 


Medium 


Lai'ge 


NIST 


1 


4s^ 


'So 


0.0000 


0.0000 


0.0000 


0.0000 


2 


4s 4p 


3Po 


0.7789 


0.7712 


0.8130 


0.8179 


3 




3P, 


0.7912 


0.7859 


0.8255 


0.8325 


4 




3P2 


0.8179 


0.8188 


0.8529 


0.8654 


5 




'Pi 


1.2842 


1.2353 


1.2567 


1.2004 


6 


4p2 


'D2 


2.0516 


1.8880 


1.9396 


1.9428 


7 




3Po 


1.8912 


1.9003 


1.9485 


1.9300 


8 




3P> 


1.9072 


1.9207 


1.9655 


1.9509 


9 




3P2 


1.9303 


1.9585 


1.9973 


1.9922 


10 


4s 4d 






2.2974 


2.3557 


2.3468 


11 








2.2996 


2.3577 


2.3488 


12 




^^D3 




2.3029 


2.3608 


2.3518 


IP 






4.8651 


4.8835 


4.9191 


5.0211 


Se'+ 


Index 


Config. 


Term 


Small 


Medium 


Large 


NIST 


1 


4s 


'Si/2 


0.0000 


0.0000 


0.0000 


0.0000 


2 


4p 


-Pl/2 


1.0065 


1.0063 


1.0272 


1.0276 


3 




-P3/2 


1.0521 


1.0520 


1.0840 


1.0795 


4 


4d 


-D3/2 




2.5111 


2.5749 


2.5773 


5 








2.5168 


2.5826 


2.5835 


I.P 




5.8878 


5.8863 


5.9592 


6.0052 



Se'' 



Index 


Config. 


Term 


Small 


Medium 


Large 


NIST 


1 


3dl» 


'So 


0.0000 


0.0000 


0.0000 


0.0000 


2 


3d''4s 


^D3 


4.0123 


4.0164 


4.0173 


4.0073 


3 






4.0351 


4.0393 


4.0401 


4.0289 


4 




'D, 


4.0731 


4.0772 


4.0781 


4.0712 


5 




'D2 


4.1138 


4.1179 


4.1187 


4.1035 


6 


3d'4p 


-^P2 


5.0873 


5.0914 


5.0929 


5.0913 


7 




-^F3 


5.1359 


5.1400 


5.1415 


5.1360 


8 




3P1 


5.1468 


5.1509 


5.1524 


5.1545 


9 






5.1720 


5.1762 


5.1777 


5.1771 


10 




'F2 


5.1760 


5.1801 


5.1816 


5.1809 


11 




3Po 


5.1813 


5.1854 


5.1869 


5.1922 


12 




'D2 


5.2463 


5.2505 


5.2520 


5.2408 


13 




'F3 


5.2465 


5.2507 


5.2522 


5.2448 


14 




^^D3 


5.2703 


5.2745 


5.2760 


5.2705 


15 




'Pi 


5.2816 


5.2858 


5.2873 


5.2696 


16 






5.3055 


5.3098 


5.3113 


5.2989 


17 






5.3128 


5.3170 


5.3185 


5.3104 


I.pl" 












11.4209 



" The internal identification scheme within AUTOSTRUCTURE likely mislabeled the Se^^ 4s^ 4p4d and 4s 4p^ 'Dt terms, leading to the large 
discrepancies in calculated vs. experimental energies. If the identifications are exchanged, the energy agreement improves significantly. 
Since we did not compute the electronic structure of Se^*, it was not possible to calculate the ionization potential of Se*^. 
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Table 4. Calculated Einstein A-coefficients (in s"') compared to literature values. A number in parentheses indicates a power of ten (i.e., 5.0(- 
10) = 5.0x10-'"). 



Se" 


i 


j 


Small 


Medium 


Large 


Previous^ 


1 


2 


1.50(-1) 


1.64(-1) 


1.64(-1) 


1.71(-1) 


1 


3 


1.53(-4) 


4.53(-4) 


4.29(-4) 


1.96(-4) 


2 


3 


9.28(-3) 


7.38(-3) 


8.58(-3) 


1.03(-3) 


1 


4 


7.59(-l) 


7.84(-l) 


7.97(-l) 


6.40(-l) 


2 


4 


1.50(-1) 


1.47(-1) 


1.51(-1) 


1.08(-1) 


3 


4 


5.71 (-4) 


6.25(-4) 


5.28(-4) 


7.86(-5) 


1 


5 


4.50(-2) 


1.26(-2) 


2.46(-2) 


1.89(-1) 


2 


5 


6.83 


6.88 


7.26 


7.61 


4 


5 


5.26(-l) 


5.47(-l) 


8.44(-l) 


2.44 


1 


6 




7.41(5) 


7.57(5) 


1.70(6) 


2 


6 




1.36(5) 


1.44(5) 


3.24(5) 


1 


7 




7.92(8) 


7.92(9) 


2.13(8) 


2 


7 




4.00(8) 


4.01(8) 


9.79(7) 


3 


7 




1.30(8) 


1.31(8) 


3.30(7) 


1 


19 




4.42(5) 


1.26(7) 


1.32(8) 


1 


20 




9.38(6) 


1.30(8) 


1.46(7) 


1 


21 




1.69(8) 


6.61(8) 


5.28(8) 


Se+ 


(' 




Small 


Medium 


Large 


Previous" 


1 


2 


3.21(-1) 


2.72(-l) 


1.97(-1) 


3.50(-l) 


1 


3 


2.33(-2) 


2.32(-2) 


1.84(-2) 


1.88(-2) 


2 


3 


2.13(-3) 


1.66(-3) 


1.10(-3) 


2.78(-3) 


1 


4 


1.61 


1.63 


1.24 


2.30 


2 


4 


5.87(-l) 


7.10(-1) 


5.31(-1) 


1.23 


3 


4 


2.30(-2) 


5.09(-2) 


4.00(-2) 


9.37(-2) 


1 


5 


3.54 


3.64 


2.78 


4.95 


2 


5 


1.10 


1.24 


9.13(-1) 


2.24 


3 


5 


5.67(-l) 


6.98(-l) 


5.22(-l) 


1.25 


4 


5 


2.67(-3) 


2.32(-3) 


1.57(-3) 


4.91(-3) 


1 


9 




1.37(9) 


1.52(9) 


9.96(8) 


1 


10 




1.41(9) 


1.56(9) 


1.03(9) 


1 


11 




1.50(9) 


1.65(9) 


1.14(9) 


Se^+ 


(' 


J 


Small 


Medium 


Large 


Previous" 


1 


2 


5.72(-2) 


6.75(-2) 


9.19(-2) 


9.59(-2) 


1 


3 


9.63(-5) 


1.15(-4) 


1.84(-4) 


1.66(-4) 


2 


3 


1.03(-1) 


1.23(-1) 


1.48(-1) 


1.35(-1) 


1 


4 


1.81 (-3) 


1.36(-3) 


2.27(-3) 


2.37(-4) 


2 


4 


6.21(-1) 


7.34(-l) 


8.70(-l) 


8.00(-l) 


3 


4 


1.15 


1.32 


1.50 


1.21 


2 


5 


1.22(1) 


1.57(1) 


1.69(1) 


1.70(1) 


3 


5 


1.42(-1) 


3.38(-l) 


3.62(-l) 


5.81(-1) 


4 


5 


7.46(-l) 


2.17 


1.57 


2.81 




i 


J 


Small 


Medium 


Large 


Previous" 


1 


2 


5.51(-1) 


8.06(-l) 


7.57(-l) 


7.58(-l) 


1 


6 


2.78(8) 


2.39(8) 


2.66(8) 


2.79(8) 


2 


6 


2.06(7) 


1.14(7) 


1.33(7) 


4.90(7) 


2 


7 


2.62(8) 


2.14(8) 


2.41(8) 


2.98(8) 


Se*+ 


(' 


7 


Small 


Medium 


Large 


Previous" 


1 


3 


2.30(6) 


4.14(6) 


4.13(6) 


1.24(7) 


1 


5 


8.04(9) 


7.23(9) 


8.50(9) 


1.79(10) 




(■ 


j 


Small 


Medium 


Large 


Previous" 


1 


2 


2.31(9) 


2.31(9) 


2.56(9) 


2.22(9) 


1 


3 


2.64(9) 


2.64(9) 


3.01(9) 


5.07(9) 


2 


4 




9.67(9) 


1.05(10) 


1.95(10) 


3 


4 




1.76(9) 


1.88(9) 


1.80(9) 


3 


5 




1.07(10) 


1.14(10) 


1.64(10) 



" References for forbidden tr a nsitions i n the ground configuration: Se" ^Biemont & Hansen" 1986a"), Se+ and Se-+ fBiemont & Hansen" 1986b|), 
Se^^ (iBie mont & Hansen|[l987h.lMortonl (l2000i) is the reference for Se°-Se^^ pe rmitted transitions, while those for Se*^ and Se'^ are taken from 
iMigd alek & Stanek' jl989l) iBahr et alj jl982h . and lCurtis & Theodosioij ( 119891) . To our knowledge, tabulated Se*"* Einstein A-coefficients have 
not been published. 



